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Multivariate central limit theorems for 
averages of fractional Volterra processes and 
applications to parameter estimation 

Ivan Nourdin* David Nualartt and Rola ZintouF 


Abstract. The purpose of this paper is to establish the multivariate normal 
convergence for the average of certain Volterra processes constructed from a 
fractional Brownian motion with Hurst parameter H > p Some applications 
to parameter estimation are then discussed. 


1 Introduction 

Let B H be a fractional Brownian motion with Hurst parameter H >i In 
this paper, we deal with fractional Volterra processes i = 1,..., k, of the 
form 


Xi(t) 



s)dB H (s ), 


t > 0, 


( 1 . 1 ) 


where ay : [0, oo) —> M are measurable functions satisfying suitable integra- 
bility conditions (to be precised later on). 

The special case of k — 1 and x\(u) = ae~ eu , with a,6 > 0, corresponds 
to the fractional Ornstein-Uhlenbeck process, which may be defined as the 
unique solution to the stochastic differential equation 


f X(t) =-QX(t) + aB H {t), t> 0 
\ X(0) = 0, 


In (11.21) and everywhere else, the dots over X and B H are used to indicate 
differentiation with respect to t. More generally, consider the following p-th 
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(1.3) 


order stochastic differential equation driven by B H : 

f XW(t) =Ej=l\e,X (j \t) + aB H {t), t> 0 
\ A(0) =... = X(P“ 1 )(0) = 0. 

In (11.31) . the superscript ^ denotes j-fold differentiation with respect to t. 
Then, it can be proved that the p-dimensional process (X, X^\ ... ,X <J> ~ 1 " > ) 
is of the form (11.11) with k = p for suitable functions Xi, i = 1,..., k. 

Considering that data are sampled from an underlying continuous-time 
process, the solution X to (ll.3|) can, for instance, serve as a model for (possi¬ 
bly irregularly spaced) discrete time long memory data. In such a situation, 
parameters 9 0 ,..., 6 k _i are usually unknown and, therefore, they must be 
accurately calibrated from the observation of X. This is how we naturally 
arrive to the issue of showing a central limit theorem (CLT) for the parameter 
estimators in the model (11.31) . in order, e.g., to construct confidence intervals. 
For instance, in the situation where those estimators are obtained by means 
of the method of moments, we may be naturally led to show a multivariate 
CLT for random vectors taking the form 

j-L^ T [(X^(t)P -E(x^\t)) m ‘}dt-. j = 1., 

for given powers mi,... ,rrik G N* and differentiation indices G 

{0,... ,p — 1}. We refer to [UBifiJ for different types of central limit theorems 
for parameter estimators in this type of models. 

Motivated by these statistical problems, our purpose is to derive general 
central limit theorems for functionals of the process X solution to (11.11) . More 
precisely, let : R —> M, i = 1,..., k, be real mesurable functions satisfying 

f fi(x)e~ x2 / 2 dx = 0 and f f 2 (x)e~ x2 ^ 2 dx < oo. (1.4) 

Jr Jr 

The second condition in (ll.4|) ensures that /) can be expanded in Hermite 
polynomials, namely 

OO 

fi = ^ a i,iHi with < °°, (!-5) 

1=0 

whereas from the first one we deduce that a^o = 0. The first goal of the 
present paper is to answer the following question. 

Question A: As T —y oo, can we exhibit reasonable conditions ensur¬ 
ing that a multivariate CLT holds for the family of random vectors Ut = 


2 










(U h T,...,U kiT )? Here 


U,,T = -j= £ f, {^) it , with = JWWl, (1.6) 

and X t . i — 1,..., k are the fractional Volterra processes solution to (11.11) . 

Since, in general, the processes are not stationary, we stress that one 
cannot directly apply the classical Breuer-Major theorem (see Theorem 12.dp 
to positively answer Question A. Nevertheless, following the approach devel¬ 
oped in Nourdin, Peccati and Podolskij (5] (see also [U Chapter 7]) we prove 
the following result. 

Theorem 1.1. Let q, denote the Hermite rank of fi, that is, the smallest 
value of l such that the coefficient a^i of Hi in HI., 51) is different from zero. 
Set g* = mini^fc g* and assume that g* ^ 2. Consider Ut = (C/pr, • • •, £4 ,t), 
where U i: T is given by H1.6\) . If H e (|, 1 —^-) and if the functions x, defining 
X t satisfy both 

/ r 

I / \xi(u)xj(v)\ |v — u — a\ 2H ~ 2 dudv ] da < oo (1.7) 
\J[ 0,oo) 2 / 



and 


r]i a H(2H — 1) I Xi(u)xi(v)\v — u\ 2H ~ 2 dudv G (0, oo), (1.8) 

V >/[0,oo) 2 


for all i, j = 1,..., k, then 


U T l ™N k ( 0,A) asT^ oo, 

where A = (Ay)^^ is given by 

A _ V' i\dI l (2H — l) 1 

/ , a i,l a j,V- 11 


l=qiVqj 


vin 


i 


X 


(1.9) 


( 1 . 10 ) 


Xj(n)xj(v)|v — u — a\~ H 2 dudv ) da 


[0,oo) 2 


That one must divide by a quantity depending on t in (11.61) . namely cq(i), 
may appear to be not very convenient for applications. This is why we also 
address the following related problem. 
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Question B: Can one find constants f t > 0, as well as suitable assump¬ 
tions on fi, Xi and H , so that Vr = (V\,t, ■ ■ ■, 14, t ) satisfies a CLT? Here 


Kt = 



Xj(t) 

6 


dt. 


( 1 . 11 ) 


Whatever the value of fi, observe that the variance of X t (t )/Q is different 
from 1 for most of the values of t. For this reason and because Hermite 
polynomials are the orthogonal polynomials associated with the standard 
Gaussian distribution, it seems difficult to deal with general functions fi 
while trying to answer Question B. This is why we restrict our analysis to 
the situation where /* are polynomials , which is not a loss of generality for 
the applications we have in mind (see Section [4]). More precisely, we have 
the following result, which provides a positive answer to Question B. 

Theorem 1.2. Suppose that fi — Pi, i — 1,..., k, are real polynomials and 
denote by g* the Hermite rank of Pi. Set g* = min^j^gj and assume that 
g* ^ 2. Consider Vr = (Vi,t,..., 14,t) given by \1.11\) . where fi = ry is given 
by (E3). If H E (|, 1 — 2 ^-) and if the functions Xi defining Xi satisfy ( | 1.7\) , 
STB as well as 

I \xi(u)xi(v)\ ((u A v) V l) \v — u\ 2H ~ 2 dudv < oo, (1-12) 

J [0,oo) 2 

then 

V T '2? N k ( 0, A) asT -> oo, (1.13) 

with A still given by M.lOi) . 

In the last section of our paper, we discuss an application of Theorem 
P to the problem of parameter estimation in the fractional CAR (k) model, 
which generalizes the model introduced in |3j. 


The paper is organized as follows. Section 2 contains preliminary results 
and concepts. The proof of the two main results, namely Theorems 11.11 and 
11.21 is then provided in Section 3. Finally, an application of Theorem 11.21 is 
discussed in Section 4. 


2 Preliminaries 

2.1 Fractional Brownian motion 

Throughout the paper, B H = (B ff (t )) fg r denotes a fractional Brownian mo¬ 
tion (fBm in short) with Hurst index H E (|,1), defined on a complete 


4 














probability space (fkJ 7 , P). That is, B H is a zero mean Gaussian process 
with covariance 

E[B«B«] = i (|t| 2 " + | S | 2 * - 1 1 - sH . 

We further assume that the a-field T is the completion of the cr-held gener¬ 
ated by B H . We denote by the closure of the space of step functions on R 
endowed with the inner product 

(1m, Wa = E[(B» - B«)(B« - B«)\, 

for any a < b and c < d. We know that the Hilbert space T) is isometric to 
the Gaussian space spanned by B H , and we denote this isometry equivalently 
by x —> B H (x) or by x i —y Jx(s)dB H (s). 


2.2 Wiener integral against fBm 

ft is well-known (see, for instance, Nualart P2 Chapter 5]) that any measur¬ 
able function x : R —?• R satisfying 


\x(u)x(v)\\v — u\ 2H 2 dudv < oc 


(2.14) 


belongs to the space £), that is, it can be integrated with respect to B H . In 
this case, the Wiener integral B H (x ) = j^^x^dB 11 (s) satisfies the follow¬ 
ing isometry property: 


B H (x)\ — \\x\\^ — H(2H — 1) / x(u)x(v)\v — u 


2H-2 


dudv. 


Notice that condition (12.141) can be equivalently rewritten as 


(|x| * |x|)(t)|t|‘ i 2 dt < oo, 


(2.15) 


with x(t ) = x(—t) and where x*y denotes the convolution of two nonnegative 
or integrable functions x, y : R —?• R: 

(x*y)(t)= / x(u — v)y(v)dv. 


The following estimate will be needed in the proof of Theorem 11.11 
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Lemma 2.1. For every x,y : [0, oo) —y R + satisfying the condition ( 2. If ) 
(extended with x(u) = y{u) = 0 whenever u < 0) and for any a e M, we have 

\ 2 

x(u)y(v)\v — u — a\ 2H ~ 2 dudv 


< 


' [0,oo) 2 

/ x(u)x(v)\v — u\ 2H ~ 2 dudv 

'[0,oo) 2 


y{u)y(v)\v — u\ 2H 2 dudv. 


' [0,oo) 2 


Proof We have, by Cauchy-Schwarz inequality, 


du x(u) / dv y(v)\v — u — a 


2H-2 


dux(u ) / dvy(v + a)\v — u H 2 

1 




H 2 (2H — l) 2 
1 


E 


x{u)dB H {u ) / r/(w + a)dB H (v) 


H 2 {2H~iy E 

/ x(u)x(n)|n — u\ 2H ~ 2 dudv 
J [0,oo) 2 


/ r°° \ 2 


ft ' 00 \ 2 

( / x{u)dB H (u) ) 

E 

( / y{v + a)dB H (v) ) 

o / 


\J-a J 


y(u)y(v)\v — u\ 2H 2 dudv. 


[0,oo) 2 


□ 


2.3 Hermite polynomials and Wiener chaoses 

For any integer 1, we denote by $)® p and fj 0p , respectively, the pth tensor 
product and the pth symmetric tensor product of f). The pth Wiener chaos 
of X, denoted by PL p , is the closed linear subspace of L 2 (Q) generated by 
the random variables {H p (B H (x)), x G Sj, ||x||^ = 1}, where H p is the pth 
Hermite polynomial defined by 

The mapping I p fx® p ) = H p (B H (x)) provides a linear isometry between lj 0p 
(equipped with the modified norm \/pl|HL® P ) an d (equipped with the 
L 2 (Ft) norm). 


2.4 Fourth moment theorem 

The following result, known as the fourth moment theorem, is a combination 
of the seminal results of Nualart and Peccati [7] and Peccati and Tudor 
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[8]. Given a sequence of random vectors in a fixed Wiener chaos whose 
covariance matrix converges, the fourth moment theorem provides necessary 
and sufficient conditions for the convergence to a normal distribution. We 
refer to Nourdin and Peccati [4j for an extensive discussion on this theorem, 
including quantitative versions obtained by means of Stein’s method and a 
wide range of applications and developments. 

In the fourth moment theorem, a crucial role is played by the notion of 
contractions. Let {e k , k ^ 1} be a complete orthonormal system in fj. Given 
/ G f) &p , g G and r G {0,..., p A q}, the rth contraction of / and g is 
the element of $)®G+ q ~' 2r ) defined by 

OO 

f ®r g = (fi e h®---® e ir)s)®r®(g,e h ®...®e ir )si®r. ( 2 . 16 ) 

r L\ y — X 

Theorem 2.2 (Fourth Moment Theorem). Let k ^ 2 and q k ^ ^ qi ^ 1 

be some fixed integers, and consider a family of kernels 

{(/l,T, • • • , fk,T )}^>0 

such that fj T G -fj 09i for every T > 0 and every j — 1,..., k. Assume further 
that 

V^ E [l qi (fi,T)I qj (fj,T)\ = A ij, V 1 ^ i, j < k. 

Then the following two conditions are equivalent: 

(i) For every i = 1,..., k and every p — 1,..., q^ — 1, 

lim || f itT (g) p /i,r||i>®a( ei -p) = 0; 

T—> oo 

(ii) as T —$■ oo, the vector (I qi (fi T ),I gk (f kT )) converges in distribution 
to the k-dimentional Gaussian vector N k ( 0, A). 

2.5 Breuer-Major theorem 

We conclude this preliminary section with a continuous-time version of the 
celebrated Breuer-Major CLT for stationary Gaussian sequences. 

Theorem 2.3 (Breuer-Major). Let (X(t))t^o be a zero mean stationary 
Gaussian process with unit variance, and let f : M —» M be a measurable 
function satisfying f R f 2 (x)e~ x ! 2 dx < oo. Let us expand f in terms of Her- 
mite polynomials, namely 

OO 

/ = aiHi with o < °°- 
1=0 
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Suppose that = 0 and let q denote the Hermite rank of f (that is, q is the 
smallest value of l such that the coefficient ai of Hi is different from zero). 
Finally, assume that J R \p(t)\ q dt < oo, with p : M —> M the autocovariance 
function associated with X, that is, 

E[X(s)X(t)] = p(t — s), t,s^ 0. 

Then, as T —» oo, 



(2.17) 


where a 2 = fm p(t) l dt G (0,oo). 


3 Proofs of the main results 


3.1 Proof of Theorem 11.11 


For the sake of clarity, the proof of Theorem 11.11 is divided into several steps. 

Step 1: Going away from zero. We claim that, in order to prove (11.91) it 
is enough to show that, for at least one fixed To > 0, one has 


law 


U T ™N k ( 0,A) as T —)• oo, 


(3.18) 


where Ut = (U^t, • • •, f4,r), with 



(The only difference between Uj.t and Ujp is that the integral defining the 
former is between 0 and T.) Indeed, using among other the Hermite ex¬ 


pansion (II. 5p of each fj and then Cauchy-Schwarz inequality to obtain that 
\E[X j (t)X j (s)]\ ^ aj(t)aj(s), we get 



E[Xfit)Xfis)] 



as T —> oo. Hence, if (13.181) holds true, then (EL9j) takes place as well. 



















Step 2: Checking that the covariance matrix of Ut converges whenever fj 
are Hermite polynomials. We can write, for any j and as t —> oo, 


<Jj(ty = H(2H — 1) / Xj(u)xj(v)\v — u\ 2H 2 dudv 




—> H(2H — 1 ) / Xj{u)xj{y)\v — u\ 2H 2 dudv = p 2 . 

J [0,oo) 2 

Let us choose T 0 > 0 large enough so that <Jj(t) 2 ^ |r/| for all t ^ T 0 and all 
j G {1,..., k}. We shall check the convergence of the covariance matrix of 
Ut with this To at hand and when fj = H Pj , j — 1,..., k, for given integers 
Pj ^ qj. In this case, one has, for any T > T 0 and any i, j G {1,..., k}, 


Pi 


mru^T] = n r „ r ,)y 


^-{pi=Pj} p I db 


T 

fb—To 


'Po,T] : 


fE[x,( s )xyt)]y dsdt 
V (Tiisjajit) J 


da 


T 0 Jb-T 


{ E[Xj(b)Xj(b — a)] 
V (Ti(b)aj(b- a) 


Pi 


P', 


i r T 


fb—To 


^-{Pi=Pj} J1 ' db 


da 


x 


JT 0 Jb-T 

H(2H - 1) 


ai(b)aj(b - a) 7[ 0 ,&]x[o,6—a] 
= l {pi=Pj} Pi\ / <L(T,a)da, 


Xj(u)Tj(n)|n — u — af 2 dudv 


Pi 


(3.19) 


where 


<h(T, a) — l{| a |s:T-To} 


/»TA(a+T) 


'ToV(a+To) 


LT(2hr-i) r 

(Ti(b)(Tj(b - a) J [ 0) 6] x [ 0 , 6 —£ 


Xj(M)xj(n)|n — u — a\ 2 2 dudv\ db. 

J 


A 


Pi 


Using first that <7j(t) 2 ^ yrj 2 for all t ^ T 0 and all i (by dehnition of T 0 ) and 
then Lemma [2.11 we deduce that, for any a G M. and any T > T 0 , 


1 fe=Pj}l^( T ’ a )l 

/ 2H(2H — 1) 

V ViVj 
f H(2H — 1) r 

\ ViVj J[ 0,oo) 2 


< 2 P 


|xj('u)xj('u)| l w — u — a\ 2H 2 dudv 


' [0,oo) 2 


A 


Pi 




\Xi{u)Xj{v) \ \v -u- a 


)qiV qj 


(3.20) 
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Moreover, due to the fact that <Ti(t) —y rj l > 0, one has, as T —> oo 

r lA(f + l) 




r t q g+T Q 
T v T 


7[0,6T]x[0,feT-a] - U - a| 2g ^ 




ai(bT)(jj(bT — a 
(H{2H - 1) f 


' [0,oo) 2 


Xi(u)xj(v)\v — u — a\ 2 ~dudv 


\ ViVj 

By dominated convergence (see also (1 1. 71) ). one deduces that 

F ,p P 1 . 1 

H Ui,TUj, T \ -> 1 tK-n) - ~pHm - 

'll 'Ij 


\Pi 

) 1 • 




(3.21) 


X 


A 


Pi 


/ Xi(u)xj(v)\v — u — af z dudv I da 

'[0,oo) 2 / 


Observe that (I3.2ip coincides with after taking into account that _/) = H Pi 
and fj = H Pj . 


Step 3: Proving \3.18\) whenever fj are Hermite polynomials. To do so, 
we shall make use of the Fourth Moment Theorem 12.21 As in the previous 
step, suppose that fj = H p ., j = 1,..., k, for some pj ^ qj. One then has 
Uj, T = I Pj ( 9j,t ), with 


Sj,T 


i r T P ® Pj 

—j= [ — yy—dt. 

VT JTo 


(3.22) 


In (I3.22lh ej )t is a short-hand notation for the function u i —> Xj(t — u)l[o ) t](u). 
According to Theorem 12.21 to conclude that (13.181) takes place we are left to 
check that, for all j G {1 ,..., k} and all r E {1,... ,Pj — 1}, 

\\gj,T djgW —t 0 as T — y oo. (3.23) 


Let us compute q h T <SV 9j,r■ We hnd 


9j,T 9j,T 


1 

T 



e j,tl 


o ®{pj~r) 

' J ,*2 




ElXjttjXjfaWdhdti. 


As a result, using moreover that <jj{t) 2 ^ \rf- for T > T 0 and introducing 
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z,m = j‘ , | Xj[t — s)\dB f (after extending B H to the whole M), we obtain 
II 9j,T ®r gj,T II 2 

= A f , , ,, / 77' ,, mx^x^YEix^x^uw 

T J [To,T] 4 ( cr j(^l) cr j(^2)c r j(^3)c r j(^4)) Pj 

XE[X j (t 1 )X j (t 3 )]^- r E[X j (t 2 )X j (tY)] Pj ~ r dt 1 dt2dt 3 dU 

4 ?j r 

EIZ^Z^YEIZ^Z^Y 


< 


rj^Po T 2 jy 0iT ]4 


x J E?[Z i (t 1 )Z J -(i 3 )]w- r J E?[Z i (t 2 )Z i (t4)]^- r dtidi 2 di3di4- 

(3.24) 


Set t| = i?[Zj(0) 2 ] and let hj t T denote the function 


h i> T ~ j>i 


r/VT, 


ef y t J dt, where e j>t (u) = \xj(t - 


Expressing the right-hand side of (13.241) by means of || hj t T <8V %r|| 2 leads to 

(3.25) 


4 Pi r 4pj 

II 9j,T Z>r 9j,r\\ 2 A -Til—I I'.i-T ®r ^j,t|| 2 - 


4 Pj 

dj 


On the other hand, it is straightforward to check that Zj is a stationary 
Gaussian process and that 

From assumption (II.7|) and with Lemma 2.1, we deduce that f R \pZj(t)\ P:i dt < 
oo. Indeed, 


PZj (t) — H{2H — 1 ) / \xj(u)xj(v)\\u — v — t\ 2H 2 dudv. 

J [0,oo) 2 

As a consequence, Breuer-Major Theorem 12.31 implies that I Vj [hj^ ) con¬ 
verges in law to a Gaussian. According to Theorem 12.21 one deduces that 
|| hj )T <g) r hj tT \\ —> 0 as T —> oo, implying in turn that (13.231) holds true (see 
(I3.25p h and thus completing the proof of (11.91) in the particular case where 

fj = H Pr 


Step 4: Proving /II..91) whenever fj are polynomials. More precisely, let us 
suppose in this step that, for each j = 1 ,,k one has 


fj 'y ' Uj.jHi, 

i= qj 


(3.26) 
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for a finite integer rrij. Owing to the Cramer-Wold device, it is actually im¬ 
mediate to apply the conclusion of Step 3 in order to get (11.9[) in the case 
where fj given by (13.261) . 


Step 5: Proving /II..91) in all generality. Finally, let us consider the general 
situation of fj given by (11.51) . To reach the conclusion in this case, and taken 
into account Step 4, it remains to show that, for any fixed j = 1,..., k, 



Using identity (13. 19[) and its associated bound (13.20)) . one has, for any fixed 
j — 1,... ,k: 


E 

1 

T 


/ 1 oo 

yrk 


V(t) 

a,(t) 


dt 


’[To,T]^ l= 


EH 


a jKMs) 


dsdt 


< — 




g[V(t)V( s )] 

T J[T 0 ,T] 2 \ a j{t) a j( s ) 

oo 

J2 Ua h 


qj oo 

dsdt x lla^ i 


l=m 


X 


(2H(2H - 1) 

V Vi 


[0,oo) 2 


\xj(u)xj(v)\\v — u — a\ 2 ~dudv I da, 


<h 


which, thanks to (11.71) and Lemma 2.1, tends to zero as m —?• oo. 


□ 


3.2 Proof of Theorem 11.2 


Set 


Yj tt = j Xj(t — s)dB H (s) and Lj,t= f Xj{t — s)dB H {s), 
J— oo J — OO 


so that Yj = Xj + Lj. It is straightforward to check that Yj is stationary and 
that rfij = Var(y j (0)). 

We have rjj = Var (Yj(t)) for all t by stationarity. So, since P has Hermite 
rank qj and since (11.71) takes place, by applying the same arguments than 
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the ones used in the proof of Theorem 1 1.1 1 but with Yj instead of Xj, we can 
prove that, as T —> oo, 


1 

vf 






™N k ( 0,A). 


Thus, to reach the desired conclusion it suffices to show that, under (\\.\2\ 
and for any fixed j = 1,... ,k, and for any integer p ^ 1, 


1 

VT 



as T —» oo. 


(3.27) 


We shall divide the proof of (13.27|1 into two steps. 


Step 1: We shall show that E L“} k t dt is bounded for any integer k ^ 
1. Indeed, one can write, with p 2 k denoting the k-th even moment of the 
standard Gaussian, 


E / L%dt 


= p2kH k (2E[ — l) k f dt ( f dudvXj(t — u)Xj(t — v)\v — u 
do \J (— oo.Ol 2 


2H-2 


< P2kH k (2H-l) h 

< P2kH k (2H-l) k 


dt 


' [t,oo) 2 


dudv \xj(u)xj(v)\\v — u 


2H-2 


' [0,oo) 2 


dui . . . dU 2 k \Xj{ui)\ . . . I Xj{u 2 k)\ 


X \U2 — U\ 


\2H-2 


U2k ~ u 2 k-i\ 2H 2 x mini'll, i = 1 ,..., 2k}. 


Now, using that minjwj, i = 1,...,2 k} ^ Yli=i iihn{n 2 i-i, u 2 i}, and taking 
into account condition (11,12)) . we deduce 


E / Lf t dt 


< kp 2 kH k (2H-iy 


[0,oo) 2 


dudv \xj(u)xj{y)\\v — u 


2H-2 


k-1 


X 


' [0,oo) 2 


dudvu A v \xj(u)xj(v)\\v — u 


2H-2 


< kp 2 kH k (2H-iy 


[0,oo) 2 


dudv ((w A v) V l) \xj(u)xj(v)\\v — u 


2H-2 


oo. 
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Step 2. Let us observe that 


1 

vf 




(3.28) 


For any fixed p ^ 1 and any 1 ^ k ^ p, one has, using the Cauchy-Schwarz 
inequality, Step 1 and that Y is stationary, 



where 0 < p < 1. So, by letting T —> oo and then p —> 0, the desired 
conclusion (j3.27j) follows, thus concluding the proof of Theorem 11.21 □ 


4 An application to the estimation of param¬ 
eters in the fractional CAR(/c) model 


Consider the fractional CAR (k) model, that is, the solution X to: 


fc-i 

X {k \t) = J20iXW(t) + aB H (t), t > 0. (4.29) 

i=0 

Here, X^ indicates the kth derivative of the solution process A" and 0j are 
real parameters considered as being unknown. Moreover, up to appropriate 
scaling it is not a loss of generality to assume that cr = 1. 

In the sequel, we are going to illustrate the use of our Theorem 1 1.21 to the 
estimation problem in (j4.29j) . To keep the things as simple as possible, we 
shall only consider the case k = 2, we shall assume zero initial conditions for 
X and we shall put some restrictions on 0 O , 0i. 

More precisely, let X be defined as the unique solution to 

X(t) = e 0 X(t) + 01 X(t) + B H (t), AT(0) = AT(0) = 0, (4.30) 


with 9 q , 6 i < 0 and 9\ + 40 o > 0. The roots of the characteristic equation 
r 2 — 9pr — 9 0 = 0 are 


V = 


0i + \/ 9\ + 40o 
2 


and 


0i — \/0i + 40o 
2 


(4.31) 
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The solution X to (I4.30P is given by 


p P p(t-s) _ P q(t-s) 

X(t) = / ---- dB H (s ). 

Jo p-q 

The processes X\ = X and X 2 = X are of the form (11.11) . with the corre¬ 
sponding functions 


Xi(t) = 


e Pt _ e <?* 

p-q 


and X 2 (t) = x\(t) = 


pe pt — qe qt 
p-q 


We shall apply Theorem 1 1.21 to X\ = X and X 2 = X and to the polynomials 
P\(x) = P 2 (x) = x 2 — 1, which have Hermite rank 1. Since p and q are 
negative numbers, the functions x\ and x 2 satisfy conditions (11.71) (11.81) and 
(1 1 .12 p . As a consequence, Theorem 11.21 implies the following convergence in 
law as T tends to infinity, provided H £ (|, |): 


Vf ^ T (X(f) 2 , X(t) 2 )dt - moo j 4 N 2 {0, A), 


(4.32) 


where rrioo = (r]f,pl), with r/,;, i — 1,2, defined in (11.81) . and A is the covari¬ 
ance matrix appearing in (11.101) . In order to explicitly compute rrioo and A, 
we shall use a Fourier transform approach. 

Computation of moo-' The first component of the vector m 00 is given by 


pl = H(2H -l) / (ai *£i)(t)|t| 


2H-2 


dt. 


The Fourier transform of X\ is given by 
= — (^-^7 

p-q \p-K q~K 

On the other hand, the Fourier transform of \t\ 2H ~ 2 is K 2 H- 2 \f,\ l - 2H , for some 
constant K 2 H- 2 - Therefore, using Plancherel theorem we can write 


Vi = 


d 


H 


(p-q) 2 

d-H 

(p-q) 2 

e H 

(p-q) 2 


p-if q-if 
1 1 


Kl 


1-2 H 


+ 


df 

2 


\e~ 2H df 


p 2 + £ 2 q 2 + £ 2 pq + i 2 , 

e H (\p\~ H -\q\~ H ? 


{\p\~ 2H + \q\~ 2H - 2(pq)~ H ) = 


(p - qy 


15 

































where e H = 2d H f 0 °° and d H = k 2H _ 2 H(2H - 1). 

For the second component of we can write 


Vl = E 


f ‘°° pe pt — qe qt 


dB H (t) 


= H(2H — 1) / (x 2 *x 2 )(t)\t\ 2H 2 dt. 


d P-Q 
The Fourier transform of x 2 is given by 


^2(0 = 


P 


Q 


p-q\p-i£ Q-% 


Therefore, 


dl = 


d 


H 


(p - q) 2 
dH 

(p - q) 2 

e H 

(p - q) 2 


P 


p-i£ q-if 


P 


+ 


lei l ~ 2H di 

2 pq 


p 2 + f 2 q 2 + f 2 pq + C 


{\p\ 2 - 2H + \q\ 2 - 2H -2{ Pq y- H ) = 


1*1 l ~ 2H dt 


e H (\p\ l - H -\q\ l - H Y 

C v-qf 


Computation of the matrix A: From (11.101) . taking into account that Pi = 
P 2 = H 2 , we obtain 


= 2H 2 (2H— iy / (xi*Xj)(u)(xi*Xj)(v)\u + a\ 2 2 |u + a 


1 2H—2 


dudvda. 


We know that 


\u + a 


2H-2 


v + a 


2H-2 


da = ku\u — v 


4H—3 


for some constant kn■ Therefore, using again the Fourier transform, we can 
write 


A ij = k H 2H 2 (2H — l) 2 / (xi * Xj)(u)(xi * Xj)(v)\u — v 


4H-3 


dudv 


= a H / \Txm 2 \Fx 3 {f)m- iH dy 

J R 

where an = ku2H 2 (2H — l) 2 KiH- 3 - The previous computations lead to the 
following expression for the components of the matrix A. 

Ai,i = t-WiMpI- 1 -" + lor 1 " 4 " + 


(p - q) A 


+(3h 


an 

2 


p 2 — q 2 


(kl 


1-4 H 


\P\ 


1-4 H\ 


p z — pq 


-((pq) 


h~ 2H 


\P\ 


1-4 H\ 


q z — pq 


-((pq) 


t-2 H ill—4 H 
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A.,2 = r 4 1 i(a»(|pr- 4 ' , + l«r-“ + 4(p 5 )J- M ) 


and 


(p-q) 
+(3h 


p2+92 (kr-“-i P r-“)-2d±«(w 


p2 _ g2 


5-2 H 


p- — pq 


\p\ 


1-4 H\ 


q z — pq J J 


^2.2 ~ 


( P - q y 

w 2fq 


OiH 

2„2 


3-4 H 


+ \q\- 4H + *(pq)-^ M ) 


-4-2# \ 


(kr-“ - w 1 - 1 ") - ipq 


p2 _ g2 


p z — pq 


((pq)->~ 2M ~ \ P \ ) 


Apq A 

q 2 — pq 

In the above expressions the constants an and (3h are given by 


-^-((™A 2H - ki 1 - 4 "))). 

q z — VQ / / 


OL}i — 2 CLh 


r°° f 2 ~ AH df 

L (i+eF 


p0 ° i 2 ~ AH di 

i + e ■ 


and (3h = 2a 


# 
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